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Observation of B Meson Decays to rjn and rj K 



Abstract 

We present preliminary measurements of the B meson decays B + — > t]tt + , B + — > rjK + , and B° — ► 
r]K . The data were recorded with the BABAR detector at PEP II and correspond to 88.9 x 10 6 BB 
pairs produced in e + e~ annihilation through the T(4S) resonance. We find the branching fractions 
B(B+ -» T/7T+) = (4.2±J;jj ± 0.3) x lO" 6 and B(B + -» r,K + ) = (2.8±g;f ± 0.2) x 10~ 6 . We set a 
90% CL upper limit of B(B° —* r/X°) < 4.6 x 10~ 6 . The time-integrated charge asymmetries are 
A ch (B + 777T+) = — 0.51±g;fg ± 0.01 and Aa(# + -> rjK+) = -0.32^;^ ± 0.01. 
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1 Introduction 



We report the results of searches for B decays to the charmless final states 1 f?7r + , rjK + , and r]K°. We 
reconstruct the 77 mesons in both of the dominant final states 77 — > 77 (r/ 77 ) and 77 — ► 7r + 7r~7r (773^). 
The if is reconstructed as ifg — > 7r + 7r~. For the charged modes we also measure the direct CP- 
violating time-integrated charge asymmetry, A c h = (r - — r + )/(r~ + where T ± = T(B ± — > 
77/^). 

The interest in these decays was sparked by the first reports of the observation of the decay 
B — > 77' if [1] in 1997. It had been pointed out by Lipkin six years earlier [2] that interference between 
two penguin diagrams and the known 77/77' mixing angle conspire to greatly enhance B — > 77' if and 
suppress B — ► 77X . Due to a parity flip for the vector if*, the situation is reversed for the B — > 77'if * 
and -B — > 7/if * decays. Though the general features of this picture have already been borne out by 
previous measurements and limits, the details and possible contributions of singlet diagrams will 
only be tested with the measurement of the branching fraction of all four (77, 77') (K, K *) decays. 

It was pointed out more than 20 years ago (before the discovery of the B meson) by Bander, 
Soni and Silverman [3] that penguin loop diagrams allow substantial ( > 20%) charge asymmetries 
in some B decays, and an example they gave was B — > r]K. The necessary ingredients are to have 
two interfering diagrams with different weak and strong phases. More recently, it was pointed out 
that such charge asymmetries can be enhanced in B — > 77'^ and B — ► 77^ where the decay rate is 
small but penguin-tree or penguin-penguin interference is possible [4, 5]. A series of quantitative 
predictions have been made in the past decade with various factorization approaches [6, 7, 8, 9]. 
There is general agreement that modes such as B — > r/K, B — > rjir, and B — > 77'^ are expected to 
have charge asymmetries of 20% or larger. Most, but not all, of the quantitative calculations predict 
negative values for all three decays. A recent paper [10] shows that branching fraction and charge 
asymmetry measurements for B — ► 77'^ and B — ► r]n allow for the determination of the strong phase 
difference between tree and penguin amplitudes and the CKM angle alpha. 

The current knowledge of these decays comes from published measurements from CLEO [11] 
and conference results from BABAR [12] and Belle [13]. Table 1 summarizes these previous results. 
We present here analyses incorporating new data. 

2 Detector and Data 

The results presented in this paper are based on data collected in 1999-2002 with the BABAR 
detector [14] at the PEP-II asymmetric e + e _ collider [15] located at the Stanford Linear Accelerator 
Center. An integrated luminosity of 81.9 fb -1 , corresponding to 88.9 million BB pairs, was recorded 
at the T{AS) resonance ( "on-resonance" , center-of-mass energy \fs = 10.58 GeV). An additional 
9.6 fb -1 were taken about 40 MeV below this energy ("off-resonance") for the study of continuum 
backgrounds in which a light or charm quark pair is produced instead of an T(4S). 

The asymmetric beam configuration in the laboratory frame provides a boost of (3j = 0.56 to 
the T(4S). Charged particles are detected and their momenta measured by the combination of a 
silicon vertex tracker (SVT), consisting of five layers of double-sided detectors, and a 40-layer central 
drift chamber, both operating in the 1.5-T magnetic field of a solenoid. Photons and electrons are 
detected by a CsI(Tl) electromagnetic calorimeter (EMC). 

Charged-particle identification (PID) is provided by the average energy loss (dE/dx) in the 
tracking devices and by an internally reflecting ring-imaging Cherenkov detector (DIRC) covering 

1 Except as noted explicitly, we use a particle name to denote either member of a charge conjugate pair. 
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Table 1: Summary of branching fraction results for B decays to r] mesons from CLEO [11], previous 
BABAR [12] measurements, Belle [13], and the present analysis. The results for all fits are given as 
well as a 90% CL upper limit if the measured yield is not judged to be significant. The overall yields 
and efficiencies (e) are given as the sum of yields and efficiencies from the two rj decay channels. 
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the central region. 

3 Event Selection 

Monte Carlo (MC) simulations [16] of the signal decay modes and of continuum and BB back- 
grounds are used to establish the event selection criteria. The selection is designed to achieve high 
efficiency and retain sidebands sufficient to characterize the background for subsequent fitting. 
Photons must have energy exceeding a threshold dependent on the combinatoric background of 
the specific mode: E 1 > 30 MeV for the two photons used to reconstruct the 7r° in rj — > tt + tt^tt° 
candidates, and E 1 > 100 MeV for rj — ► 77. Additionally, we require that the cosine of the center 
of mass decay angle for r/ 77 daughters, relative to the flight direction of the 77, have an absolute 
value of less than 0.95 (0.97) for neutral (charged) decays involving the r/ 77 . 

We select ?? 77 , i] 37T , and 7r° candidates with the following requirements on the invariant mass in 
MeV/c 2 of their final states: 490 < < 600 for r/ 77 , 520 < m n < 570 for 7737,-, and 120 < m n o < 
150. For Kg — ► 7T + 7r _ candidates we require 488 < rriK s < 508. Typical resolutions are about 
16 MeV/ c 2 for ?? 77 , 4.5 MeV/ c 2 for r/ 37r , 2.8 MeV/ c 2 for K° s , and 7MeV/ c 2 for tt°. For K° s candidates 
we require that the three-dimensional flight distance from the event primary vertex be > 2 mm, 
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and the two-dimensional angle between flight and momentum vectors be < 40 mrad. 

We make several particle identification (PID) requirements to ensure the identity of the signal 
pions and kaons. Tracks in 773^ candidates must have DIRC, dE/dx, and EMC responses consistent 
with pions. For the charged B + — ► r]K + decay, the prompt charged track must have an associated 
DIRC Cherenkov angle between —5 a and +2 a from the expected value for a kaon. For B + — > r/7r + , 
the DIRC Cherenkov angle must be between —2a and +5 a from the expected value for a pion. 

A B meson candidate is characterized kinematically by the energy-substituted mass toes = 



y (^s + Po -pb) 2 /Eq — p 2 B and energy difference AE = E* B — \yfs, where the subscripts and 

B refer to the initial T(4S) and to the B candidate, respectively, and the asterisk denotes the 
T(4S) frame. The resolutions on these quantities measured for signal events are 30 MeV and 

3.0 MeV/c 2 , respectively. We require \AE\ < 0.2 GeV and 5.2 < toes < 5.29 GeV/c 2 (the lower 
limit is 5.22 GeV/c 2 for r/ 77 vr+). 

3.1 Tau, QED, and continuum background 

To discriminate against tau-pair and two-photon background, we require in 773^ channels that the 
event contain at least five (four) charged tracks for neutral (charged) B pairs. In 7/ 77 analyses, we 
require three (two) tracks for neutral (charged) B pairs. 

To reject continuum background, we make use of the angle Qt between the thrust axis of the 
B candidate and that of the rest of the tracks and neutral clusters in the event, calculated in 
the center-of-mass frame. The distribution of cos 9t is sharply peaked near ±1 for combinations 
drawn from jet-like qq pairs and is nearly uniform for the isotropic B meson decays; we require 
I cos^tI < 0.9. A second B candidate satisfying the selection criteria is found in about 10-20% of 
the events. In this case the "best" combination is chosen as the one closest to the nominal r/ mass. 

The remaining continuum background dominates the samples and is modeled from sideband 
data for the maximum likelihood fits described in Section 4. 

3.2 BB background 

We use Monte Carlo simulations of B°B° and B + B~ pair production and decay to look for possible 
BB backgrounds. Most BB backgrounds in these analyses come from other charmless decays. 
From these studies we find no evidence for significant BB background in the 77 — > 7r + 7r~7r° decay 
chains. 

For the 77 — > 77 modes we find potential BB backgrounds from several charmless final states, 
which we treat with additional event selection criteria. To reduce background from 7r°7r + , ir°K + , 
and 7t K°, we eliminate r; 77 candidates that share a photon with any 7r° candidate having mo- 
mentum between 1.9 and 3.1 GeV/c in the T(4S) frame. Additionally, we remove high-energy 
photons to suppress background from K*j, by requiring E 1 < 2.4 GeV. We find a small remaining 
BB background in rj^ir (tj^K) from r\p (rjK*). To discriminate between these and the signal we 
include a BB component in the likelihood fits for modes with 77 — > 77, as described in Section 4.1. 

4 Maximum Likelihood Fit 

We use an unbinned, multivariate maximum-likelihood fit to extract signal yields for our modes. 
A sample of events to fit is selected as described in Section 3. 
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4.1 Likelihood Function 



The likelihood function incorporates four uncorrelated variables. We describe the B decay kine- 
matics with two variables: AE and ?7?es- We also include m v and a Fisher discriminant T which 
describes energy flow in the event. The Fisher discriminant combines four variables: the angles 
with respect to the beam axis, in the T(4S) frame, of the B momentum and B thrust axis, and the 
zeroth and second angular moments Lo,2 of the energy flow about the B thrust axis. The moments 
are defined by 

Lj = ^Pi x |cos0i| J ' , (1) 

i 

where 6i is the angle with respect to the B thrust axis of track or neutral cluster i, pi is its 
momentum, and the sum excludes the B candidate. 

As measured correlations among the observables in the selected data are small, we take the 
probability distribution function (PDF) for each event to be a product of the PDFs for the separate 
observables. We define hypotheses j, where j can be signal, continuum background, or (for modes 
with ?7 77 ) BB background. The product PDF (to be evaluated with the observable set for event 
i) is then given by 

V) = VMes) • Vj(AE) • Vj{Tl ■ Vj(m v ). (2) 
The likelihood function for each decay mode is 

exp (- Yi) 

N i - IIE*^. ( 3 ) 

* 3 

where Yj is the yield of events of hypothesis j found by the fitter, and ./V is the number of events 
in the sample. The first factor takes into account the Poisson fluctuations in the total number of 
events. 



4.2 Signal and Background Parameterization 

We determine the PDFs for signal and BB background from MC distributions in each observable. 
For the continuum background we establish the functional forms and initial parameter values of 
the PDFs with data from sidebands in toes or AE. We allow several background parameters to 
float in the final fit. 

The distributions in m v , and in toes an d AE for signal, are parameterized as Gaussian func- 
tions, with a second or third Gaussian as required for good fits to these samples. Slowly varying 
distributions (combinatoric background under the rj mass and AE peaks) are parameterized by 
linear functions. The combinatoric background in toes is described by a phase-space-motivated 
empirical function [17]. We model the T distribution using a Gaussian function with different 
widths above and below the mean, and include a linear contribution of 1-3% in area to account 
for outlying events. The linear term ensures that the significance of the signal is not overestimated 
relative to background. Because of the rarity of outlying events this component is not particularly 
well determined in some data samples, but we have checked that the yield and its significance are 
insensitive to choices of a linear component over the conservative range 1-6%. 

We check the simulation on which we rely for signal PDFs by comparing with large data control 
samples. For toes and AE we use the decays B~ — > ir~D°, D° — > K~ir + ir°, which have similar 
topology to the modes under study. For we use inclusive resonance production. 
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5 Fit Results 



By generating (from PDF shapes) and fitting simulated samples of signal and background, we 
verify that our fitting procedure is functioning properly. We find that the minimum ln£ value in 
the on-resonance sample lies well within the ln£ distribution from these simulated samples. 

The efficiency is obtained from the fraction of signal MC events passing the selection, adjusted 
for any bias in the likelihood fit. This bias is determined from fits to simulated samples, each equal 
in size to the data and containing a known number of signal MC events combined with events 
generated from the background PDFs. We find biases ranging from 1% to 4%, depending on the 
mode. 

Table 2: Final fit results for B + — > rjh + and B — > r]K° , where r\ — > 7r + 7r~7r° and r\ — > 77. We 
report branching fractions for the two 77 decay channels separately (£>) and after combining the 
results of the two channels (Combined B). Systematic contributions are included in the signficance 
values. The Corrected B for the charged modes is the branching fraction after correcting for 
crossfeed from one charged mode into the other. 



Fit quantity 






r/ 77 7r + 






rj 11 K + 






Fit sample size 


















On-resonance 


9477 




6933 


5383 




5884 


1270 


1435 


Off-resonance 


1104 




1168 


630 




959 


158 


183 


Signal yield 
On-res data 
Off-res data 

Selection e (%) 


28.0+i° 8 ° 
1 1+ 2 - 1 

1 - 1 -1.2 

23.3 




39.611J;? 
28.7 


0.6±I;| 
22.6 




34.31U 
30.6 


2.6+fi 

o +0 - 7 

22.6 


o-o±8:g 

24.8 


UBi (%) 


22.6 




39.4 


22.6 




39.4 


7.8 


13.5 


exUBi (%) 


5.2 




11.3 


5.1 




12.1 


1.76 


3.34 


Stat. sign, (cr) 


4.3 




5.7 


2.4 




5.7 


0.8 


3.2 


B(xW- 6 ) 
Combined B 
Stat. sign, (cr) 


e.oifi 

4.5 


J - y -1.0 

+1.0 , 

-0.9 ^ U - J 

7.0 


O + 1.8 

J - z -1.5 

3.2 


9+O.9 
°- z -0.8 

t° 8 7 ±0.2 
6.2 


1 7 +2.6 9 q+1.6 
l.<_ 2 .o z - y -1.3 

2.6lg;| ± 0.2 
3.3 


Corrected B 


4 9 +1 -° 


±0.3 


2-818:? 


±0.2 


2 6 +0 - 9 
z -°-0.8 


±0.2 



90% C.L. UL(incl. syst.) < 4.6 



BkgA ch -0.00 ±0.01 -0.02 ±0.01 -0.02 ±0.01 -0.00 ± 0.01 

Signal Ah -0.50 ±0.31 -0.51 ± 0.24 -0.56 ±0.55 -0.25 ± 0.26 
Combined A ch — 0.51±g;?g ± 0.01 — 0.32±g;?| ± 0.01 

Stat. sign, (cr) 2.5 1.4 



In Table 2 we show the results of the fits for off- and on-resonance data. Shown for each 
decay mode are the number of events that were fit, the signal yield, the efficiency (e) and daughter 
branching fraction product (J[ Bi), and the central value of the branching fraction. We also show the 
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branching fraction results after combining the two r/ decay channels, before and after a correction 
for crossfeed between the two charged channels (see Section 7), and the statistical significance of 
this combined result. For r]K° we quote a 90% CL upper limit. The statistical error on the number 
of events is taken as the change in the central value when the quantity — 21n£ changes by one 
unit. The statistical significance is taken as the square root of the difference between the value 
of — 21n£ for zero signal and the value at its minimum. For the charged modes we also give the 
charge asymmetry A c h- 

In Fig. 1 we show projections of m^s an d AE made by selecting events with signal likelihood 
(computed without the variable shown in the figure) exceeding a mode-dependent threshold that 
optimizes the expected sensitivity. 




m ES (GeV/c ) 



0.1 0.2 

AE (GeV) 



Figure 1: Projections of the B candidate tties an d AE for B + — * r]n + (a, b), B + — > r/K + (c, d), 
and B° -> rjK° (e, f). Points with errors represent data, shaded histograms the r/ 77 subset, solid 
curves the full fit functions, and dashed curves the background functions. These plots are made 
with a cut on the signal likelihood and thus do not show all events in the data samples. 
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6 Systematic Uncertainties 



Most of the systematic errors on yields that arise from uncertainties in the values of the PDF 
parameters have already been incorporated into the overall statistical error, because their back- 
ground parameters are free in the fit. We determine the sensitivity to parameters of the signal PDF 
components by varying these within their uncertainties. The results are shown in the first row of 
Table 3. This is the only systematic error on the fit yield; the other systematics apply to either the 
efficiency or the number of BB's. 

The uncertainty in our knowledge of the efficiency is found to be 0.8iVi%, 2.5iV 7 %, and 3% 
for a Kg decay, where N t and iV 7 are the number of signal tracks and photons, respectively. 
We estimate the uncertainty in the number of produced BB pairs to be 1.1%. The estimate of 
systematic bias from the fitter itself (1-2%) comes from fits of simulated samples with varying 
background populations. Published world averages [18] provide the B daughter branching fraction 
uncertainties. We account for systematic effects in cos#t (1%) and in the PID requirement (0.5%) 
on the prompt charged track. Values for each of these contributions are given in Table 3. 

A study of the charge asymmetry as a function of momentum for all tracks in hadronic events 
bounds the tracking efficiency component of charge-asymmetry bias to be less than 1%. Samples of 
B and D*-tagged D — > Kir decays provide additional crosschecks that the bias is small. We assign 
a systematic uncertainty for A c h of 1.1% based on the tracking study and a small PID contribution 
determined from the D* studies. 

Table 3: Estimates of systematic errors (in percent) for the B + — > rjh + and B° — ► r]K° modes. We 
specify which systematics are uncorrelated (U) or correlated (C) between rj decay channels. 



Quantity 




7/ 77 7T + 




ri~ f ~ f K + 


V3*K° 




Fit yield (U) 


3.9 


3.7 


8.4 


4.5 


20.7 


2.3 


Fit efficiency /bias (U) 


1.9 


1.3 


1.3 


0.8 


1.0 


1.7 


Track multiplicity (C) 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


Tracking eff/qual (C) 


2.4 


0.8 


2.4 


0.8 


3.7 


2.1 


7r%/7 eff (C) 


5.0 


5.0 


5.0 


5.0 


5.0 


5.0 


K° s efficiency (C) 










2.9 


2.9 


Number BB (C) 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


Branching fractions (U) 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


MC statistics (U) 


1.1 


1.1 


1.0 


1.1 


1.1 


1.0 


cos# T (C) 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


PID (C) 


1.4 


1.0 


1.4 


1.0 


1.0 




Total 


7.5 


6.9 


10.5 


7.3 


22.0 


7.2 


Uncorrelated 


4.6 


4.2 


8.6 


4.8 


20.8 


3.2 


Correlated 


6.0 


5.5 


6.0 


5.5 


7.2 


6.4 
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7 Combined Results 



We next combine the r\ — > 3ir and rj — > 77 branching fraction measurements. We do this by 
first forming for each rj decay mode the convolution of C from the fit with a Gaussian function 
representing the uncorrelated systematic error. The curves — 21n£ are shown in Fig. 2, for each 
r) mode and for their sum. For the time integrated charge asymmetries the corresponding — 21n£ 
plots are given in Fig. 3. 

The results at this stage are given in the row labeled "Combined £>" in Table 2. For the charged 
modes we must apply a correction for kaon-pion crossfeed arising from imperfect PID. In studies 
with kaon and pion samples tagged kinematically from the decays D* + — > ir + D°, D° — > K~ir + we 
find that 9 ±2% of pions are accepted by our kaon selection and vice versa. After correcting for this 
and adding the associated systematic uncertainty we obtain the final measurements summarized in 
Section 8. 

8 Conclusion 

We report preliminary measurements of branching fractions and A c h for B meson decays to rj with 
a charged kaon or pion, as well as the branching fraction for B° — > 7]K°. We find statistically 
significant signals in the charged B decays. The branching fractions are 

B(B + 7/7T+) = (4.2±£;g ± 0.3) x 10" 6 , 
B(B + -► V K + ) = (2.8±g;| ± 0.2) x 10~ 6 . 

For the neutral B decay, we find B(B° -> rjK°) = (2.6l°;l ± 0.2) x 10~ 6 . Since the statistical 
significance of this result is only 3.3a, we determine a 90% CL upper limit: 

B(B° -► V K°) < 4.6 x 10~ 6 . 

These results supersede the previous BABAR measurements [12]. Our measurements of the CP- 
violating charge asymmetries in the charged modes are 

Ach(B+ -+ 7/7T+) = -0.51±g;?§ ±0.01, 
A ch (B+^r)K+) = -0.32±g;?| ±0.01. 

These charge asymmetry results are in agreement with the theoretical expectations discussed in 
Section 1 and rule out substantial positive asymmetries. 
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B(B + ^r|Jt + ) B(B + ^r|K + ) 
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2.5 5 7.5 10 



B(B°^riK ) 

Figure 2: Distributions of — 21n£ vs branching fraction for rjir + , 7]K + and 7]K° decays. Two 
secondary channels (dashed and dotted lines) are combined to produce a final result (solid line). 
The dashed line corresponds to 77 — > 7r + 7r~7r° decays, while the dotted line corresponds to rj — > 77. 
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Figure 3: Distributions of —2 ln£ vs A c k for ??7r + and r]K + decays. Two secondary channels (dashed 
and dotted lines) are combined to produce a final result (solid line). Dashed line corresponds to 
rj — > 7T + 7T~7r decays, while the dotted line corresponds to rj — > 77. 
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